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Many theoretical studies were devoted in the past to ion-specific eflects, trying to interpret a 
large body of experimental evidence, such as surface tension at air/water interfaces and force mea- 
surements between charged objects. Although several mechanisms were suggested to explain the 
results, such as dispersion forces and specific surface-ion interactions, we would like to suggest an- 
other source of ion-specificity, originating from the local variations of the dielectric constant due to 
the presence of ions in the solution. We present a mean-field model to account for the heterogeneity 
of the dielectric constant caused by the ions. In particular, for ions that decrease the dielectric 
constant we find a depletion of ions from the vicinity of charged surfaces. For a two-plate system, 
the same effect leads to an increase of the pressure in between two surfaces. Our results suggest that 
the effect of ions on the local dielectric constant should be taken into account when interpreting 
experiments that address ion-specific effects. 



I. INTRODUCTION 

The first interest in ion specific effects dates back to 
the end of the 19th century, when Franz Hofmeister and 
his co-workers [l| measured the thermodynamic prop- 
erties of protein precipitation in various salt solutions. 
They tested numerous salt species and found that the 
protein solubihty properties can be arranged in a specific 
ionic order. This classification of cations and anions is 
known today as the Hofmeister series. Since then many 
other experiments demonstrated that ion-specific effects 
take place in a large variety of chemical and biological 
systems 0. These experiments include measurement of 
surface tension at air/water interface proteins sta- 

bility @] , forces between charged surfaces such as mica or 
silica surfaces 0, II] , osmotic pressure in systems contain- 
ing biological macromolecules [Qhllj and many more. 

Noticeable efforts have been devoted over the years to 
understand the physical mechanisms that lead to ion- 
specific effects 0, [12h15{ . A typical problem is to ob- 
tain density profiles of ions near charged surfaces such 
as macromolecules, membranes and colloids or neutral 
dielectric interfaces {e.g., air/water interface), and use 
the profiles to calculate macroscopic quantities such as 
osmotic pressure and surface tension. 

Theoretical endeavors offer a variety of perspectives 
regarding the origin of the ion-specific effects. Most re- 
cently, dispersion interactions [l^, [3 depending on the 
polarizability of the ions, have been proposed to furnish 
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the missing link between ionic profiles and ion-specific 
interactions. These interactions add an additional term 
to the total ionic energy close to a non-charged dielec- 
tric interface. This term varies as B/z^, where z is the 
distance of the ion from the surface, located at z = 0. 
The ion-specificity then emerges from a variation in the 
specific ion polarizability, leading to different values of 
the coefficient B. 

Apart from the dispersion interactions, ion-hydration 
interaction has been suggested to lead to an effec- 

tive short-range ion-specific interaction with the bound- 
ing surface. For example, an attractive potential of the 
ions to the surface in the shape of a square well, takes into 
account their larger affinity as compared to ion-water in- 
teraction. This attraction then boosts the concentration 
of the ions in the surface vicinity. 

In a similar fashion but with an opposite effect, On- 
sager and Samaras dealt in their seminal work [l7] with 
repulsive image charge interactions as a source of ionic 
depletion from a non-charged dielectric interface. Image 
charge interaction is a short-range surface interaction, 
having an exp(— 22;/Ad)/z dependence, where Ad is the 
Debye screening length. This repulsive interaction be- 
comes significant for air / water interfaces and is responsi- 
ble for the increase of surface tension of saline solutions 
as compared with pure water. More recently, additional 
interactions have been proposed to augment the original 
Onsager and Samaras model, and offer an explanation 
to ion-specific surface tension of electrolyte solutions, in 
agreement with experiments [H, [l^ . 

In another approach Q it was suggested that surface 
charge density is a crucial parameter controlling ionic 
specificity. The density is determined by the net charge 
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and specific volume (i.e., radii) of tfie ions. For small 
ions, the surface charge density can be large, leading 
to a strong attractive interaction with the polar water 
molecules. This results in a "hard" hydration shell with 
high energy cost of removing it. Oppositely, for large 
ions, the surface charge density is smaller, and there- 
fore, the energy cost of breaking the ion- water complex is 
lower, leading to a "soft" hydration shell. Due to the dif- 
ferent energy costs of "hard" and "soft" hydration shells, 
it is favorable to have pairing where both cations and 
anions are either small or large. An interesting feature 
of this ionic pairing mechanism is a possible explanation 
of a reversal in the order of the Hofmeister series, as ob- 
served in several experiments 0- 

Besides the analytical approaches, computer simula- 
tions gain much needed insight into the intricacies of the 
ionic specificity (2^, [2l| , and are used to investigate var- 
ious cross-interactions between the system components 
(solutes, solvent and surfaces). Computer simulations 
capture a very detailed microscopic picture as compared 
to continuum theories such as the Poisson-Boltzmann 
(PB) theory and its extensions. However, they often lack 
the simplicity and predictive power of analytical models 
and their straightforward intuition. 

From the current state of theoretical and experimen- 
tal studies, it appears to be quite clear that ion-specific 
effects are due to the interplay between the ion-ion, ion- 
solvent and ion-surface interactions, taking into consid- 
eration both electrostatic and non-electrostatic interac- 
tions. However, it is still not well understood what are 
the important physical and chemical properties of the 
solutes and solvent that determine these specific interac- 
tions. 

The traditional PB theory accounts for electrostatic 
interactions on a mean-field level with well understood 
drawbacks and limitations [22| . It neglects all ion-specific 
effects except for the ion valency, which means that it 
treats identically all ions of the same valency. Apart from 
the electrostatic interactions, other non-electrostatic in- 
teractions can be added into the theory. This was done in 
the standard DLVO theory of colloid stability , where 
the total interaction decomposes into a sum of the van 
der Waals and electrostatic interactions. 

In the present work we would like to suggest a phe- 
nomenological approach to treat different ionic species. 
Our approach shares some of the simplicity of the PB 
theory, while taking into account the ionic specificity on a 
mean-field level. In particular, we focus on the effect the 
ions have on the local solution polarity as codified by its 
static dielectric function. The static dielectric function is 
treated as a spatially-dependent function s{z), where the 
spatial dependence is given implicitly by the local ionic 
concentration, thus introducing ionic specificity in an im- 
plicit manner. In the next section we elaborate on the 
dependence of the dielectric constant on the ionic concen- 
tration and discuss its physical range. In section iHll we 
present a mean-field model that accounts for the effect of 
the ions on the dielectric constants. Finally, results for 



TABLE I: Values of the linear coefficient /3, Eq. dTJ, in units of 
M^^ for various salts, [adopted from the references as indicated]. 
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density profiles and the inter-plate pressure are presented 
in sections HVl and IVl respectively. 



II. ION DEPENDENT DIELECTRIC 
CONSTANT 

The static dielectric function of an electrolyte solution 
is generally found to be smaller than that of the pure 
solvent [24, 25] . This decrement of the dielectric constant 
has been attributed to various sources, which underly the 
changes in the dielectric response of the solution. We 
mention here two of the most important ones: 

Ionic polarizability. Each ion in an aqueous solu- 
tion creates a cavity and displaces one or several water 
molecules in this process. Since the ions usually have 
smaller static and dynamic polarizabilities than water, 
they modify the dielectric response of the solution [2(i |. 
This effect, though omnipresent, is nevertheless small and 
does not have substantial consequences. 

Hydration shell. The gist of the dielectric decrement 
in aqueous electrolytes is connected with the structural 
modification of water molecules in the immediate prox- 
imity to the ion due to large electrostatic fields emanating 
from the dissolved ions 26] . By the action of these strong 
electrostatic fields, hydration shells are formed around 
solvated ions, where nearby water molecules are oriented 
along the ion electrostatic field, leading to an additional 
pronounced dielectric decrement as is depicted schemat- 
ically in Fig. [1] 

The lowering of the dielectric constant in the vicinity 
of an ion can be attributed to the effects described above, 
which represent the quintessential mechanisms of the di- 
electric decrement. In an ionic solution the same effect 
would lead to an overall reduction of the total solution 
dielectric constant as a function of the ionic concentra- 
tion. In addition to these primary mechanisms, there are 
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Hydration Shell 



FIG. 1: A schematic drawing of water molecules around a mono- 
valent cation in presence of an electrostatic field, E. Arrows rep- 
resent the dipoles of water molecules while the hydrogen bond 
network is not shown explicitly. The water molecules in the vicin- 
ity of the cation form a hydration shell and are oriented along the 
field lines generated by the cation, leading to a decrease of their 
contribution to the screening of the external electrostatic field. 
The more remote molecules are less affected by the cation and 
orient themselves along the external field lines. Furthermore, the 
cation polarizability results in a different contribution to the elec- 
trostatic screening as compared to the contribution of a water 
molecule. 



other possible sources of the dielectric decrement such as 
excluded volume [13] and dynamical effects [l^ that can 
have a noticeable contribution. 

In numerous experimental studies [29l - l32j and simula- 
tions [ill , the dependence of the dielectric constant was 
found to depend linearly on the salt concentration n for 
molar concentrations ranging between zero and 1.5 M: 



e{n) = Ef) + I3n , 



(1) 



where e(n) is the dielectric constant of the ionic solu- 
tion, eo is the dielectric constant of pure solvent, and (3 
is a phenomenological coefficient (in units of M~^) of the 
linear term. The most common case where /3 < corre- 
sponds to a dielectric decrement. Throughout this paper 
the value of Eq — 80 is taken to be that of pure water. 

The dielectric decrement parameter /? has salt-specific 
values [25,, 29,ri33] . In Table 1 we show a range of such 
values adopted from different experiments and computer 
simulations. The cited values are all negative and vary 
from about -8 to -21 (dielectric decrement). There is also 
a noticeable spread of the /3 values, especially when one 
compares experiments with computer simulations. Ex- 
amining the /? values given in Refs. [25| and [1^ for ho- 
mologous series of the halides: CI, F and I, and the alka- 
line metals: Li, Na, K, Cs, and Rb, some remarks can be 



made. In some cases, a trend can be seen, where the mag- 
nitude of (3 decreases as the ionic radius (size) increases 
for the alkalines. This can be seen for the homologous 
series of XCl where X is the alkaline ion. For the ho- 
mologous halide series, the trend of with the halide 
ionic size is not monotonous. We note that the above lin- 
ear form of the static dielectric constant e{n) depends on 
the concentration of the salt and cannot, in general, be 
entangled in a simple manner into separate ionic contri- 
butions of the cations and anions. Any such separation 
requires additional hypotheses and experiments that are 
unavailable at present. 



III. THE MODEL 

In what follows we concentrate on the dielectric decre- 
ment effects of ions (when /? is negative), and system- 
atically assess the modification wrought in the equilib- 
rium ionic profiles as well as interactions between two 
apposed charged surfaces. Note that we treat the prob- 
lem of counter-ions only (no co-ions) in order to simplify 
the model and results. 

We set up a phenomenological model based on a mod- 
ification of the standard PB theory, which includes the 
effect of the counter-ion concentration n{r) on the lo- 
cal dielectric constant, in a consistent manner, via the 
local variation of n{r) in e{n). This model can be for- 
mulated in two distinct but equivalent versions. First, 
one can dress it in the form of the modified PB theory 
starting with the appropriate generalization of the PB 
free-energy functional where the fundamental quantity is 
the local electrostatic potential. Equivalently, as is shown 
in the Appendix, the PB free-energy functional can also 
be cast in the form of a density functional theory (DFT), 
where the fundamental quantity is the local density of 
the counter- ions, n(r). 

The mean-field free energy is similar to the regular 
PB ansatz, except for the coupling between the electro- 
static field, — VV', and the positively charged counter-ion 
concentration, n. For simplicity, we delimit ourselves to 
systems that are translationally invariant in the lateral 
(x, y) directions having a spatial dependence only on the 
z-direction. The free energy then reads: 

-f fcBTn(logn-l)-H+E.e|CT|V., (2) 

where A is the lateral cross-sectional area, e is the elec- 
tron charge, fee the Boltzmann constant, T the temper- 
ature and ^ the chemical potential. The last term is the 
contribution of the charged surfaces to the free energy, 
where — cr < is the negative surface charge density, tps 
is the electrostatic surface potential, and the last term is 
summed over all charged surfaces. Note that throughout 
the manuscript ea is the surface charge density and a is 
the corresponding number density. 
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Using the variational principle for the free energy, 
Eq. ^ , with respect to n and ip gives the Euler-Lagrange 
equations determining the electrostatic potential and 
density profiles: 



= 



Sip 



— - 
Sn 



Anen + ^ ^£(n)?/;'^ = , 

/i = ——e'{n)ip'^ + 61^ + k^T log n 
on 



(3) 
(4) 



where e'{n) = de/dn. 

As was discussed in our previous work as long 
as the system has only one-dimensional spatial inhomo- 
geneity, the pressure is a z-independent constant, and in 
fact, represents the first integral of the Euler-Lagrange 
equations. Using Eq. (25) of Ref. we arrive at the 
following form of the pressure, P 



P = -— {£{n) + e'[n)T}\ -0'^ + k^Tn = const . (5) 

The pressure is composed of an appropriately modified 
Maxwell stress tensor, which takes into account the den- 
sity dependence of the dielectric constant [11] , as well as 
the standard van't Hoff term proportional to the counter- 
ion concentration. 

Combining Eqs. ([S])-®, the following first-order differ- 
ential equation for n is derived: 



dn 
dl 



/27re2 



k^T fin) 



(6) 



the pressure in Eq. ([5]) vanishes, P 
f{n) from Eq. ([T]) is then given by: 



0. The function 



f{n) = e{n) 



e{n) + ne'{n) 



(9) 



Evaluating Eq. ([5]) at the surface with e = £s and n — rig 
we obtain an algebraic equation for n^, given by: 



27re2a2 / , de 



kBTe^in,) 



(10) 



As has already been noted in the previous section, for 
ionic concentrations of < n < 1.5 M, the function e{n) 
can be regarded to a reasonable accuracy as linear in n, 
Eq. ([T|). This is a key assumption of our model where 
all the effects discussed in Sec. |lT] contribute to a single 
phenomenological parameter /3 = de/dn. We stress that 
the linear dependence mentioned in the previous sentence 
is a phenomenological approach and cannot be justified 
theoretically for any value of concentration [stI]. 

The derivative of f{n) is given by: 



/'(") 



2/3^n' + (£o + ^n)(£o + 2/3Ti) 
2v^n(eo 4- 2/3n)3 



(11) 



and the boundary condition, Eq. (|10|) . can be expressed 
as a cubic polynomial in n^: 



I3^nl + 2/3£on' + ( £ 



ksT 



kBT 



= 0. 

(12) 



where 



A. The ^ = profile 



fin) = £(n)i 



P 



e{n) + ne'{n) 



(7) 



and P — P/k^T. The solution of Eq. ^ yields the equi- 
librium profile of the counter-ion density. The boundary 
condition at the charged interface is obtained by taking 
the variation of Eq. ([2]) with respect to the surface po- 
tential tps'- 



'd7 



— Anecr , 



(8) 



where Eg = £{ns) is the value of the dielectric constant 
extrapolated to the surface, and rig is the counter-ion 
concentration at the surface. 



IV. THE SINGLE PLATE CASE 

We first consider the case of a single charged plate 
placed at 2 = 0. This case can be regarded as a two-plate 
system in the limit of infinite inter-plate spacing, where 



The equation for the counter-ion density profile that 
we derived above, Eq. ([5]), appears very different from 
what one finds in the standard PB theory. Let us show 
that the PB solution is indeed obtained by setting (3 = 
in our model. Substituting (3 — in Eq. ([TT1) . the density 
profile obeys the following equation: 



dn 



PB 



dz 



I 87re2 
kBTeo 



3/2 
^PB ' 



(13) 



and its solution reproduces exactly the familiar profile of 
Gouy and Chapman (GC) [36| : 



1 



2TTh 



1 



Z + X: 



GC 



(14) 



where Aqc = 1/{2ttIb<j) is the GC length, Ib — 
e^/ie^kBT) is the Bjerrum length and Eg — £o for the 
standard PB case. The counter-ion concentration at the 
surface is obtained via the Graliame equation, n^^ = 
2ttIbO'^ [36| . or alternatively can be expressed in terms of 



Agc, n. 



PB 



1/(2^/bA: 



GC) 
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FIG. 2: In (a) the ion concentration at the interface, Us, is plotted as a function of /3. The inverse of Eq. (|15p is denoted by the 
solid line. The surface charge density is a = 0.01 A~^. Approximated analytical solutions for several limits are denoted by roman 
numerals (see text). In (b) the parameter q is plotted as a function of 7. The inverse of the exact relation, Eq. (|17p. is denoted by 
the solid line. The approximated analytical solutions are denoted by the same notations as in (a). 



B. The boundary condition for /3 7^ 

Returning to the general /3 7^ case, the counter-ion 
concentration adjacent to the surface, Ug, satisfies the 
more general relation given in Eq. ()12|) . Introducing the 
dimensionless parameters 7 = /3nf ^/eq and a = Us/n^^ , 
Eq. can be rewritten as: 



7 



a 



(15) 



where ^ ^ as obtained from the above mentioned 
Grahame equation. 

It is evident from Eq. p3|) that for any /3 7^ 0, 
a = Us/n^^ < 1 and the counter- ion close to the sur- 
face, is reduced with respect to PB value. Surprisingly, 
even when /3 > and the local dielectric constant in- 
creases (electrostatically favorable), there is an effective 
repulsion of the ions from the interface. 

It is also instructive to express Eq. (|15l) in terms of 
the dimensionless parameter q = Cs/eo — 1 = /Sus/so, 
which measures the /3-dependent relative change of the 
dielectric constant at the interface. The dimensionless 
algebraic equation of the boundary condition, Eq. (1151) . 
can now be written in terms of 7 and q. 



g3 + + (1 _ 2j)q - 7 - , 



and its solution yields: 



2,2 



1 + 2, 



(16) 



(17) 



In several limits the expressions of ns{p) and q{j) sim- 
plify: 

(i) For strongly negative /3 ^ — 1, equivalent to 7 <C — 1, 
the counter-ion concentration at the surface is decaying: 



'2/3 



(18) 



while the change in the surface dielectric constant ap- 
proaches a limiting value q = —0.5, leading to Ss = eo/2. 
For these values of Us and Sg the polarization energy 
[see right hand side of Eq. ([5|)] associated with the water 
molecules is equal the energy of the ions "effective" po- 
larization, yielding a limiting value for the electrostatic 
energy. 

(ii) For small \f3\ ^ 1 or I7I <^ 1, the correction to rig is 
small and given by: 



,PB 



iPB 



I3n 



PB 



£0 



(19) 



while the leading order of q is linear in 7, 5(7) ~ 7. 
(iii) In the limit of extremely polarizable ions (/3, 7 3> 1), 
Us varies as: 



' 2eon 



PB 



1 



£0 



4 V 2l3n 



PB 



and q (or e^) grows as 

9=^27(1 



4^ 



V7. 



-1/2 



(20) 



(21) 



Using Eq. (IT71) one can analyze the dependence of Ss, the 
dielectric constant value at the interface, on /?. 



One would expect that for /? > 0, where the ions increase 
the dielectric constant, the surface ionic density Us will 
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FIG. 3: Counter-ion concentration n(z) and local dielectric constant e{z) profiles for 7 — —1.8. In (a) the solid line shows the 
counter-ion concentration for negative /? = — 20M^^. The dashed line corresponds to the standard PB case (/? = 0), as given by 
Eq. I|14p. In (b) the local dielectric constant e{z) is shown for the same /? value. The surface charge density is cr = 0.01 A^^ and the 
calculated surface values are n^^ = 7.3 M, Us — 1.9 M and = 43. 



be increased as well. Surprisingly, is decreased for 
increased /3 as a result of the interplay between all the 
contributions, see [38]. 

The exact dependence of ns{/3) from Eq. (fT5|) . and 
5(7) from Eq. (|17p . is compared with the above limit- 
ing expressions in Fig. [2l There is a smooth crossover at 
712 — —1 between regimes (i) and (ii), and similarly at 
723 ~ 1 for the crossover between regimes {ii) and (Hi). 
The crossover values are taken to satisfy the conditions 
of the approximations. Calculating higher order terms 
for the three regimes gives intersection values that are 
similar to our crossover values. 

C. The profiles for /3 / 

The above analysis for the limiting regimes for 5(7) or 
equivalently Ssil), applies also for the concentration and 
dielectric profiles. We first present our results for nega- 
tive /3, which is the more prevalent case for most mono- 
valent ions (e.g., halides and alkalines) (see Table 1). 

The concentration profile is shown in Fig. [31 It has a 
plateau in the vicinity of the surface, till about distances 
of z ~ 5 A and larger, where there is an algebraic decay as 
in the regular PB profile (dashed line) . The value of 7 = 
— 1.8 used in Fig. [3] is smaller than the crossover value, 
7i2 = —1 and lies within regime (i). For completeness 
we also present in Fig. |4] the (3 > case, where 7 = 2.3 
is larger than 723 = 1 [regime (ii)]. Although the value 
of Us is reduced, the profile exhibits a regular algebraic 
decay, similarly to the PB profile (dashed line). 

In the more common case of /3 < (dielectric decre- 
ment by the ions), the ions should be depleted from 
the nearby wall region. This proximity behavior of the 



counter-ion cloud bears some resemblance to the case of 
steric effects [S^, which leads to similar (but not identi- 
cal) exclusion of counter-ions from the vicinal region. 

The crossover values 712 can be transformed into re- 
lations between /3 and the surface charge density — cr, 
yielding a regime diagram as shown in Fig. [Sfa). The 
crossover line between regimes (i) and {ii) is given by 
the relation: 

which holds only for (3 < 0. This crossover line sepa- 
rates between concentration profiles with a plateau-like 
behavior in the vicinity of the surface, to those where 
the concentration profile is nearly identical to the PB 
one. This crossover can be understood in terms of the 
dual role played by counter-ions as they accumulate at 
the surface. On one hand, just as in the standard PB 
model, the ions are attracted to the oppositely charged 
surface and diminish the local electrostatic field. On the 
other hand, for /3 < the ions reduce the local dielec- 
tric constant and, hence, are repelled from the surface. 
In regime (i), the latter prevails and leads to a plateau 
in the counter-ion density when a and |/3| are above the 
crossover line, as in Fig. [3Ka). 

Similarly, the crossover line between regimes (ii) and 
(Hi) is given by: 

which holds only for /323 > 0. In regime {Hi) large values 
of /? increases substantially the value of e close to the 
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FIG. 4: Counter-ion concentration profile n{z) and local dielectric constant e{z) for 7 = 2.3 > 723. In (a) the solid line shows 
the counter-ion concentration for positive j3 — 25 M^^. The dashed line corresponds to the standard PB case (/3 = 0), as given by 
Eq. (|14p . In (b) the corresponding local dielectric constant is shown for the same /3 value. The surface charge density is cr = 0.01 A^^ 
and the calculated surface values are = 7.3 M, ri^ = 4.7 M and £s = 197. 




FIG. 5: In (a) the regime diagram in the (/3,ct) plane in shown. Two cross-over lines separate the three limiting regimes: (/), 
(;7) and (//7). In (b) typical concentration and dielectric constant profiles for each regime are shown, with surface charge density, 
a = 0.01 A^^. The dashed line corresponds to regime (;') with /3 — — 20M~^, Us = 1.9 M and Eb = 43. The solid line corresponds 
to regime (//) with ,5 = 1 M"\ Us = 7.2 M and Es = 87. The dotted line corresponds to regime {Hi) with /3 = 25 = 4.7 M 

and Es = 197. 



surface [Fig. Hl^b)], but has no large effect on n{z) w 
npBiz). This somewhat surprising result is due to the 
attenuation of the electrostatic attraction for very large 
increase in the values of e at the surface. 

In Fig. [5jb), we demonstrate the qualitative differ- 
ence between the various regimes by plotting three cor- 
responding counter-ion concentration profiles. For all 
shown profiles, the surface charge concentration is = 
0.01 A^^, while the parameter (3 takes the values (3 = 
—20, 1 and 25 for the dashed, solid and dotted lines. 



respectively. 

For the profiles of regime (i), /3 < 0, one can calculate 
approximately the electrostatic potential in the vicinity 
of the charged surface by assuming that the counter- 
ion concentration varies slowly close to the surface and 
roughly obeys n fa Hg. The potential is then given by the 
quadratic form: 

27re 

4,iz) ~ -— {-n,z^ + 2az) . (24) 
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/ z[A] / z[A] 

FIG. 6: Profiles of electric potential, concentration and dielectric constant for parameter range as in regime (;'). In (a) the electrostatic 
potential ?/' is shown for a = 0.01 A^^ and /? — — 20M^^. The solid line corresponds to the exact numerical solution, while the 
dotted line denotes the approximated one, Eq. (|24p. and the dashed line is the regular PB (/3 = 0) result. In (b) and its inset, the 
concentration and dielectric constant profiles are shown, respectively, for the same parameters as in (a). The position of I = 2.3 A 
as calculated from Eq. I|25p is marked in (a), (b) and the inset. It corresponds to the crossover from the slowly varying behavior to 
regular PB algebraic decay. 



The width of the vicinal plateau, I, can be estimated by 
noting that the electrostatic field for which the dielec- 
tric decrement becomes negligible is simply related to 
the crossover value, ai2 given in Eq. ([22)) . By demand- 
ing that the electrostatic field at z = I matches 
— 47recri2/(eo + f^Us), the plateau width I is estimated to 
be: 



I = 2a- 



ttIbSo ' 



(25) 



where Ug ~ £o/(2|/3|) in regime (i), Eq. (ITSl) . The profile 
for z > I , can then be estimated as a regular PB profile 
with adjusted surface charge a — (T12 and a shifted z axis: 
z ^ z — I. The validity of Eqs. (p4| and ((25|) is exam- 
ined in Fig. [S] by comparing them to the exact numerical 
results. Several observations are worth noticing. First, 
the value of Z ~ 2.3 A is rather short-range. It gives an 
estimate to the width of the depleted ionic layer. Second, 
the assumption of a saturated ionic layer reproduces well 
the electric potential up to z ~ 3L 

When compared to the corresponding PB results, the 
electrostatic potential shows a marked increase depend- 
ing on the dielectric decrement. For parameters used in 
Fig. [HI the electrostatic potential almost doubles in size as 
compared with its standard PB value. The effect persists 
for z up to several dozens of Angstroms. Nevertheless, 
the more important effect of the electrostatic potential 
on the ionic density itself is much shorter ranged, and 
is attributed to the spatial variation of the dielectric re- 
sponse, see inset Fig. IHl^b). 



V. THE PRESSURE FOR THE TWO-PLATE 
SYSTEM 

We now move on to analyze the case of two apposed 
charged planar surfaces. For simplicity sake, we re- 
strict ourselves to the symmetric case where two equally 
charged plates are located at 2; = ±D/2, and evaluate 
the corresponding disjoining pressure P = P/k^T . 

For an arbitrary function e{n) and P 7^ 0, the bound- 
ary condition, Eq. ([5]), is generalized as follows: 
2^e2^2 



" kBTe^{ns) I ^ ^"^ ^ dn 



(26) 



This equation can then be cast into a form that contains 
only Us and e{ns) without its derivatives [40j . 

27reV2 , n, P 
u rr , . + log + — = const. (27) 

It shows that the pressure is in fact determined solely by 
the values of the surface energy and the (ideal) entropy 
of mixing at the surface. 

Employing the linearity assumption, e(n) = Eq + /3n, 
the derivative of f(n) is given by: 



fin) 



2p^n{n ~P) + (eo + Mi^o + ^fSn) 



(28) 



2J(n-F)(£o + 2/3n)3 



and the boundary condition, Eq. (1261) . can be expanded 
as a cubic polynomial in n^: 



/3hil + (2/3eo - /3'P)nl + {el 



2pPeo 



( 



Pel 



(29) 
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FIG. 7: The pressure P is shown as a function of the sepa- 
ration D for surface charge density a — 0.01 A^'^. The solid, 
dashed, dotted and dash-dotted lines correspond to values of 

/3 = 0, -5, -10 and -20 respectively. 



The counter- ion density profiles are given by Eqs. ([B]), 
and (Uni, while the dependence of the pressure on 
the inter-plate separation D is obtained by solving these 
equations for various values of P, and inverting the func- 
tion D{P) into PiD). 

Let us first verify that our theory has the correct PB 
limit for /3 = 0. The differential equation of the density, 
Eq. ^ reads 



dn 
dl 



' 87re2 
fceTeo 



nJ (n - Ppb) , 



(30) 



where PpB is the standard PB pressure. With K — 
K{P) = \/ 2i:l^P, the solution of the above equation 
can be shown to be 



n(z)-P 
arctauA/ ^ — —Kz, 



or can be expressed as 



n{z) = P(l + tan2 Kz) 



cos^ Kz 



(31) 



(32) 



This is exactly the standard solution of the GC equation 
between two equally charged walls. The pressure is then 
obtained by solving the boundary condition [SG^ that re- 
duces to the following transcendental equation: 



fKD\ 
KD tan f j 



D 



^GC 



(33) 



For the general case of /3 ^ 0, the basic equation for 
the counter-ion density profile as well as the boundary 



condition have to be solved numerically. In Fig.[7]the de- 
pendence of the pressure P on the separation D is shown 
for several values of /3. The dashed line of /3 — — 5M^^ 
is almost identical to the standard PB prediction, given 
by the solid hue (/3 = 0). This implies that /? = -SM'^ 
is roughly the value where the dielectric decrement be- 
gins to affect the disjoining pressure. For more nega- 
tive P = — 20M^^, the pressure strongly depends on the 
dielectric decrement and increases substantially with re- 
spect to the PB prediction. For example, at a separation 
£) = 25A, the pressure P = 0.7M with /3 = -20 M^^ is 
increased by 40% comparing to Ppb — 0.42 M. 

For small separations (and large pressure) the assump- 
tion of linear dependence of e on n breaks down due to 
high values of n > 1.5 M. The results shown in Fig. [7] in- 
clude only separations where the model is still valid. In 
the regime of extremely high surface counter-ion densities 
our model predicts unphysical results, such as negative 
local dielectric constant. This breakdown does not imply 
any deep inconsistency but is just a straightforward con- 
sequence of the linearity ansatz, Eq. This problem 
can easily be overcome by introducing a more general 
non-linear ansatz with a saturation behavior. 

The deviation AP = P — Ppb of the pressure with 
respect to the PB pressure PpB, is plotted in Fig. [SJa) 
as a function of D. The AP deviation is significant for 
separations as large as a few nanometers. For D < 30A 
the relative deviation may be as large as 65%, while for 
£) ~ 50 A the deviation is smaller than 0.05 M, leading to 
corrections smaller than 20%. In Fig.|Sl[b), the deviation 
is plotted on a log-log scale, where AP shows a power- 
law decay, AP ^ D". With a linear fit the extracted a 
exponent is a = —2.7,-2.8 and —3.1 for 13 = -5,-10 
and — 20Af~^, respectively. 

For (3-^—1 one can estimate the pressure analytically 
using rescaled separation Deft and surface charge density 
CTofi . We consider an approximated profile with a plateau 
near the surface, having a thickness: 



I 



(T - (712 



(34) 



employing the same arguments as in section HVCI Note 
that the values of Us and ai2 are now given by Eq. (|29|) 
instead of Eq. (fT2|) . The pressure associated with this 
profile can be calculated from the region where the be- 
havior is PB-like, namely, for -D/2 + I < z < D/2 - I. 
This region can be considered as an independent pro- 
file with a surface charge density Ccff = o'i2 and sepa- 
ration D(.f{ = D — 21. In analogy with the two limiting 
regimes of the standard PB theory 36], we find them 
here as well. First, ideal-gas regime where DcS <^ 
and the pressure depends inversely on the separation: 
P ~ (Jcs/Dcs, where the rescaled GC length is defined 
as = l/(27r/BO'off)- Second, when D^s ^ A^, 

the pressure dependence is given by the standard Gouy- 
Chapman result valid for sufficiently large inter-surface 
separations, i.e., P ~ 7r/(2ZB-D(,g)- 



10 



0.8 



^ 0.6 



0.4 



0.2 



J 


'. 1 1 1 


(a) 


L 
1 


■ 

\ 


- 


1 
1 

n ; 
1 


\ 
\ 
\ 

■ 




1 

J 


\ \ 




\ 


\ \ 




\ 


\ 




\ 


\ 




- \ 


\ 




\ 




\ 


'n 




\ 

1 







10 



20 



30 
D[A] 



40 



50 




10 



10 



10" 



D[A] 



FIG. 8: In (a) the pressure deviation as compared to the PB pressure (/3 — 0), AP = P — PpB, is plotted as function of D, for surface 
charge density a = 0.01 A^'^. In (b) similar results as in (a) are plotted on a log-log scale. The dashed, dotted and dash-dotted lines 
correspond to values of 13 = —5, —10 and —20 M^^, respectively. 
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FIG. 9: In (a) the total pressure P is shown as a function of the parameter j3 for surface charge density a = 0.01 A ^. In (b) the 
pressure correction, AP = P — PpB, as compared to the PB pressure, is similarly plotted. The solid, dashed, dotted and dash-dotted 
lines correspond to values of D — 40, 30, 20 and 17 A, and standard PB pressure PpB = 0.82, 0.63, 0.31 and 0.19 M respectively. 



The pressure as a function of the parameter /3 is pre- 
sented in Fig. [nja) for several values of the separation 
D. The significance of the dielectric decrement effect 
becomes substantial for (3 < — 5M~^ as can be seen in 
the plot. For smaller values of (3 the deviation from the 
standard PB value, P{(3 — 0) is small and negligible for 
all the values of D. The magnitude of the deviation 
AP as a function of /? is presented in Fig. MJ^)- For 
/3 < — 11.5M~^, the deviation is of the order of at least 
10% of the total pressure for all the values of D. For 
example, for Z? = 30A and /3 = — 11.5M^^, the pressure 
P = 0.37 M and the deviation is AP = 0.053 M. 

For completeness we present pressure profiles for /3 > 
in Fig. [TUl It is evident from comparing Fig. [5fa) and 



Fig. [To] that for positive values of l3 the deviation AP 
is much smaller. This is in agreement with the analysis 
of the density profiles (see Sec. IV. C), where a similar 
small deviation from the standard PB profile is found for 
/3 > 0. 



VI. CONCLUSIONS 

The model presented in this work accounts for local 
changes in the dielectric constant of a solution due to 
the presence of ions with an effective polarization effect. 
Assuming linear dependence of the dielectric constant e 
on the ionic concentration n, we introduce a phenomeno- 
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FIG. 10; The pressure correction, AP = P — Ppb as a function 
of D, is plotted for several /3 > values: solid, dashed and dotted 
lines correspond to /3 = 5, 10, and 20. 



logical parameter f3 that describes the relation between 
e and n: e — Eq + I3n. This parameter is ion-specific, 
and as shown in several experiments and computer sim- 
ulations, its value is negative and varies between —8 and 
-21M-1. 

The concentration dependence of the dielectric con- 
stant leads to an additional coupling between the 
counter-ion concentration and the electrostatic poten- 
tial. Consequently, the counter-ions play a dual role. 
On one hand their net charge induces an attractive in- 
teraction with the charged surface. On the other hand, 
the counter-ions effective dipole moment due to the con- 
centration dependence of the dielectric constant leads to 
another electrostatic interaction that could be either at- 
tractive for /3 > 0, or repulsive for /? < 0. The interplay 
between theses two interactions leads to changes of the 
free-energy minimal configuration. 

Analyzing the behavior in {l3,a) parameter plane, we 
find three regimes. The most pronounced effect is found 
for regime (i) and /? <C — 1 , where ions are strongly de- 
pleted from the charged surface due to the decrement of 
e and the resulting penalty of electrostatic energy. This 
depletion results in a plateau of the ionic concentration 
profile close to the charged surface. We estimate the 
plateau thickness analytically and find that it is of the 
order of few angstroms. In regime (ii) where \/3\ < 1, the 
effect of local changes of the dielectric constant is neg- 
ligible comparing to the contribution that comes from 
the net charge of the ions. The profiles in this regime 
are characterized by an algebraic decay similarly to the 
standard PB theory. For /3 ^ 1 in regime (iii), the di- 
electric constant is increased nearby the surface, leading 
to a weaker attraction of the ions to the surface. The 
surface concentration is reduced comparing to the stan- 
dard PB theory, but the overall change in the profiles is 
small. 



Calculating the disjoining pressure between two planar 
surfaces, we find a substantial increase for negative values 
of /3 of the order of —10. The deviation of the pressure 
calculated in our model as compared to the standard PB 
pressure reaches ~ 200% for /3 — —20 and a separation 
D ~ IsA. The range of this additional effect persists up 
to a few nanometers. For (3 = — 10M~^, the deviation 
is significant for separations smaller than ^ 30A. The 
dependence of the correction AP on the separation D 
exhibits an inverse power-law decay. The power varies 
between —2.7 to —3.1 for different values of /3 as discussed 
in the previous section. For positive values of /? we find 
that the pressure does not vary substantially comparing 
to the standard PB prediction. 

It is rather important to notice that the effect of the 
dielectric decrement on the disjoining pressure shows up 
as an effective "solvent structural force" [4l| , see Fig. [3 
between the two charged surfaces that acts at relatively 
small inter-plate separations, though no structural forces 
were assumed a priori in our approach. This might signal 
a more general relationship between the dielectric decre- 
ment and water structure effects. In the present formula- 
tion, where the dielectric decrement itself does not con- 
tain any spatial scale, the dependence of the additional 
"water structural force" on the inter-plate separation ap- 
pears to scale with an inverse power of this separation. 
However, for a model where the dielectric decrement lev- 
els off for larger values of the counter-ion density, the ad- 
ditional phenomenological constant will introduce a new 
length scale into the problem. Then, possibly also the 
additional "water structural force" would show such a 
length scale. 

Further refinements and applications can be consid- 
ered. For example, the model can be generalized to in- 
clude salt ions, instead of counter-ions only as was done 
here. This will introduce another screening effect (the 
Debye length) that will compete with the other contri- 
butions. Moreover, the concentration dependence of the 
dielectric constant can be taken to reproduce the real ex- 
perimental results in order to improve the quantitative 
accuracy. Finally, by generalizing the model to curved 
geometries it will be possible to calculate potential of 
mean force between spheres and cylinders. 

The results presented in this work suggest that the 
effect of ions on the local dielectric constant should be 
taken into account on equal footing as other ion-specific 
interactions such as dispersion and hydration. Further- 
more, this model may serve as a platform for a more 
detailed models, taking into account other ion-specific 
effects. 
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Appendix: The DFT formulation for two-plate 
system 

We would like to show how to extend our formula- 
tion in Sec. IIIII to the case where other interactions, be- 
sides electrostatic ones, also depend on the density pro- 
file, n{z). Just as in Sec. |Vl the equally charged plates 
are taken to be at z = ±D/2 yielding that iIj'{z = 0) = 
at the symmetric midplane, z = 0. The first Euler- 
Lagrange (EL) equation, Eq. ([3]), still remains valid 



which is exactly the second EL equation, Eq. ([4]), derived 
in Sec. IIIII within the PB formulation. 

The two formulations are indeed completely equiva- 
lent. The DFT formulation is preferable when one has 
additional terms in the free energy that depend either on 
the density or its derivatives, but this venue is left for 
future work. 



^P'iz) = 



47re 
4n{z)] JO 



n{z')dz' . 



(A.l) 



A second integration would then yield V-'(^) a function 
of n(z). Furthermore, it follows from the same mentioned 
EL equation that 



d 

dz 



47r 



47r ^ 



(A.2) 



Instead of performing the minimization on the free energy 
with respect to the density field, n, as in Eq. (jlj, we 
will first express the free energy as a functional only of 
n. To be followed only then by a free energy variation. 
Inserting Eq. (|A.2p into the square brackets of the free 
energy, Eq. ([2]), the free energy becomes 



F/A - /dz 



3^ (4?^^' 



fcs (log n — 1) — /in] + "EsCaips ■ (A. 3) 



By integrating explicitly the second term in Eq. (|A.3p it 
is easy to see that this contribution exactly cancels the 
surface term. The free energy per unit surface area then 
assumes the form 



F/A 



kBTn{logn 



1) — fin ) dz. 

(A.4) 



Substituting the expression for the local -0' field, 
Eq. (|A.1|) . we get 



F[n]/A^J 



2 e[n-M ilo^i^'W) 

kBTn{z){logn{z) — 1) — /in(z)] dz . 



(A.5) 



The free energy is now a functional of the local density 
n{z) only. The variation of the above functional with 
respect to n(z) leads to 



eip ■ 



1 de{n) 
Stt dn 



i;'^ + keTlogn- fi = 0, 



(A.6) 
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